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A 3D Cellular Automaton model developed by the authors to deal with the dynamics of N-body 
interactions has been adapted to investigate the head-on collision of two identical bound clusters 
of particles, and the ensuing process of fragmentation. The range of impact energies is chosen low 
enough, to secure that a compound bound cluster can be formed. The model is devised to simulate 
the laboratory set-up of fragmentation experiments as monitored by 4n detectors. The particles 
interact via a Lennard-Jones potential. At low impact energies the numerical experiments following 
the dynamics of the individual particles indicate a phase of energy sharing among all the particles 
of the compound cluster. Fragments of all sizes are then found to evaporate from the latter cluster. 
The cluster sizes, measured in our set-up by simulated Att detectors, conform to a power law of 
exponent ~ 2.6. 

In an attempt to duplicate the laboratory caloric curves related in particular to nuclear fragmen- 
tation processes, we introduce several temperature parameters (kinetic temperature of nucleons, 
kinetic temperature of fragments, reaction equilibrium temperatures). Theoretical caloric curves 
are then constructed for those temperature parameters we regard as physically most relevant. Our 
results show that different temperature definitions generate different curve patterns, indicating that 
the fragmentation system remains far from thermodynamic equilibrium. The pattern of the labo- 
ratory caloric curve for Au-Au collision experiments as derived from a recent analysis [NuPECC 
report Dec. 1997] is reproduced qualitatively by our reaction temperatures. 



I. INTRODUCTION 

The currently favoured theoretical approach of investigating many-particle fragmentation consists in applying sta- 
tistical mechanics to the excited compound cluster. The procedure relies on the assumption that the collection of 
particles, or compound cluster, originating from the two colliding clusters reaches a state of thermodynamic equilib- 
rium. The properties of this compound cluster can then be described thermodynamically, in terms of a small number 
of global parameters, namely the total number of particles A, possibly the total charge number Z, the total available 
energy E (regarded as measured in the centre of mass system), the total volume V occupied by the particles, and 
possibly the area of the surface E enclosing this volume (0 for a recent review of the theoretical approaches in the 
nuclear field; 0) S S S S S 111 for experiments and results on molecules and clusters of ions). In this Paper we shall 
describe a less restrictive theoretical framework applicable to a multiple fragmentation process of arbitrary nature 
(nuclear, atomic, or molecular). For our specific illustrations the relevant orders of magnitude of the global parameters 
have been chosen to capture the particular nuclear case. 

A few comments on the related laboratory experiments are therefore in order. The empirical results are collected 
essentially through the following experimental procedure. A target nucleus {At,Zt) is bombarded with a high- 
energy beam of incident nuclei {Aj, Zj) generated by an accelerator. A series of detector set-ups (such as ALADIN, 
CHIMERA, EOS, FOPI, INDRA, LASSA or MINIBALL), ideally uniformly distributed over a sphere centred on the 
target nucleus and thereby securing a Att coverage, identify the charge and mass of the coUisional fragments as well 
as their kinetic energies. 

Experiments of this nature were spearheaded in the 80s 0, 0| ■ In 0] the target was a Kr or Xe ion; the 
accelerated incident particles consisted in a beam of protons of energy in the range 80-350 MeV. The detectors 
isolated nuclear fragments in the A-range 12 < A < 31. The yield Y{A,Z) was found to be consistent with an 
expression of the form A^'^ times a Boltzmann factor depending on a temperature parameter. The exponents 
evaluated experimentally were t — 2.64 (Xe) and 2.65 (Kr). Qualitatively the outcome of these experiments was 
compared with the thermodynamic transition from a liquid phase (original target heated by the infalling beam of 
protons) to a state of formation of droplets of all sizes (fragment distribution measured by the detectors) taking 



2 



place at a critical-point temperature Tc. In fact, for the latter transition a power-law in the sizes of the droplets 
is known to hold, of exponent t = 2.33. This result is derived in the context of the mean field theory |0. Similar 
fragmentation experiments were described in |E3,^3|> which the target was an Ag, Kr, Xe or U ion; the incident 
particles were protons or carbon ions. Again fragments were detected (here in the Z-range 3 < Z < 22), which obeyed 
an approximate power law. 

The seemingly canonical interpretation of these earlier as well as of the more recent intermediate-mass ion-ion 
collision experiments (for instance the ^^^Au-^^^Au fragmentation reported on in regards the fragmentation 

process as the formal analogue of a liquid-gas phase transition occurring at a transition temperature Tc- At lower 
temperatures, T < T^, the nuclear system is a mixture of individual nucleons and fragments (liquid phase); at higher 
temperatures, T > T^, all fragments dissolve into a gas of nucleons (gas phase). 

This interpretation of an observed fragmentation process (of any nature) as a thermodynamic phase transition may 
be helpful in providing a qualitative picture for the outcome of the collision experiments. However, it should be clear 
that a many-body fragmentation process induced by violent collision cannot be regarded as a proper thermodynamic 
equilibrium phenomenon. Prior to the collision, both the target cluster, and the incident cluster, are in a stable 
internal equilibrium state. The two clusters are characterised globally by two sets of (extensive) thermodynamic 
variables {Ag, Zq,Eq,Vq, Eg), q = T,I. The energies Eg are measured in the centre-of-mass (CM) system of ions T 
and / respectively. The stability of the equilibrium of each cluster requires a zero surface pressure. At the moment 
of the incipient collision the T and / clusters merge into a single cluster C, the compound cluster C (counterpart of 
the compound nucleus C in the sense of Bohr and Wheeler jial)- The extensive thermodynamic variables are then 
essentially the sums of the extensive variables of the components T and /, Ac = At + Aj, etc, except for the difference 
in the evaluation of the energy. The energy of C, Eq, is to be measured in the CM frame of the compound system, so 
that we have Ec = Et + Ej + Eimp', the extra energy component, Eimp, is the impact energy (global kinetic energy 
of the T and / in CM of C). 

Since the impact energy is a free parameter, the surface pressure of C, the initial merger of T and /, will not 
vanish in general, and hence the compound cluster will not be in a state of thermal equilibrium. The surplus energy is 
expected to create a pressure increase (positive surface pressure), which will force C globally to expand. This expansion 
may consist in a release of clusters of particles which ultimately become independent fragments (as observed in the 
experiments), leading in turn to a contraction of a central core. Such a phenomenon of global expansion made up 
of an expanding outer envelope (together with a contracting inner core) is familiar in gravitationally bound systems 
(spherical clusters of stars; cf the classic analysis by Lynden-Bell and Wood 16]). In the general fragmentation 
problem, the central remainder C" is a new compound cluster, with fewer particles, which keeps evaporating fragments, 
thereby transforming again into a smaller central cluster C", as long as it has not reached a final equilibrium state. 
These comments make it manifest that a priori the distribution of fragments as measured by the laboratory detectors 
is not directly related to a statistical equilibrium state of a central cluster C; it is rather the system's non-equilibrium 
state that causes the evaporation of the fragments. However, provided that the central cluster evolves slowly enough, 
the component particles may have time to reach an approximate statistical equilibrium. The latter alternative is 
envisaged in gravitational systems as well. The distribution of the clusters that evaporate at a given time t is 
then expected to correspond to a chemical reaction equilibrium, or a dissociation equilibrium, at the temperature 
characterising the central cluster at time t. On the other hand, it is clear that this situation can hold only if the impact 
energy remains low enough. At higher impact energies the target will be essentially transparent to the particles of 
the incident cluster. 

The simplest statistical models by-pass the conceptual question of whether or not a statistical equilibrium holds 
in the fragmentation process. In these models the system of particles is enclosed in an energetically insulated box 
of finite volume Vb] the collection of interacting particles is then fully specified by 3 independent thermodynamic 
parameters, Vb, Aq, and Eq (total energy) which are given at the outset. Such a system is necessarily due to 
relax eventually towards a thermodynamic equilibrium. In the specific nuclear context, detailed classical molecular 
dynamics (CMD) calculations have indeed demonstrated directly the existence of a statistical equilibrium state, from 
which the thermodynamic properties of the system can be recovered Moreover, the distribution of the 

fragments with size can be evaluated in the thermal equilibrium state. Similar thermodynamic equilibrium results 
have also been derived in the framework of cellular automaton (CA) numerical experiments (in 2D [l^ and 3D 
However, the connection between these, as well as other theoretical models so far proposed [I| on the one hand, and 
the real, a priori far from equilibrium laboratory experiments on the other hand, remains entirely unclear. 

In this Paper we pursue the goal of setting up a theoretical framework capable of closely simulating the arrangement 
of the actual laboratory experiments. To this end (1) we follow in the first place the dynamics of the collision (as 
in the case of CMD calculations) without relying on an assumption of a thermodynamic equilibrium. (2) Next we 
analyse the collection of fragments which have arrived at a certain distance do from the collision site where they 
can be regarded as independent. The fragment distribution is estimated at that particular level. Our procedure 
of evaluation of the theoretical distribution thus contrasts with the conventional statistical methods which have 
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been applied by previous authors, including ourselves |l8l Il9l l20| . It conforms instead essentially to the laboratory 
readings as given by an array of detectors located at a distance do from the collision site. (3) We introduce and 
compute explicitly a variety of formal temperature parameters, namely kinetic temperatures associated with the gas 
of particles and the gas of fragments, and reaction temperatures related to the distribution of the fragments. If the 
system were in a state of genuine thermodynamic equilibrium, in particular during the initial stage of the collision 
when the two colliding clusters merge into a single compound cluster, then a true thermodynamic temperature would 
exist. Under this condition the formal temperature parameters would all be equal to the true temperature. Provided 
only that there are experimental procedures replicating the conditions of the theoretical definitions of these formal 
temperature parameters, the latter continue to provide a useful global characterisation of the system, whether or 
not a thermodynamic equilibrium state is realised: An acceptable theory must then be capable of duplicating the 
experimentally available values of these formal parameters. 

The question of the temperatures is analysed in greater detail in Section 5. The dynamical calculations (Sectiori 4) 
are carried out in the context of the Cellular Automaton (CA) model developed by the authors and discussed in [13 
and its 3D extension [23|. A major difference between the numerical experiments of [23| and the experiments of the 
present work resides in the fact that the system of component particles was confined to a finite box, while in principle 
the system of particles in our present analysis evolves in a virtually infinite lattice space. 

The detailed CA experiments we report on refer to a nuclear fragmentation involving collisions between identical 
ions (cf the ^^^Au-^^^Au collision of ALADIN, We treat these collisions as being head-on. In a future work we 

plan to extend our model to deal with collisions of non-zero impact parameter. 



II. THE CA MODEL 

The C A framework adapted to the simulation of the classical dynamics of interacting particles is discussed in |l9j | , 
for the specialised nuclear context; the particles are nucleons, and we consider low enough energies. We indicate here 
only the particular features of this model when applied to the collisional dynamics we are concerned with, and recall 
the typical orders of magnitude of the model parameters. 



A. Geometry of Lattice Space. Particle Kinematics 

Our CA universe ?7 is a cubic lattice of cells, of toroidal topology (periodic boundaries), with a typical size L = 
127. A lattice cell is identified by a position vector r of integer Cartesian coordinates r = (a;, y, z), with x, y, z taking 
the values -(L - l)/2, -(L- l)/2 + 1, . . ., -1, 0, 1, . . ., {L - l)/2- 1, (L - l)/2. An individual cubic cell has a side 
A chosen of the order of magnitude of the range of the nuclear forces (w 2 fm) . The timestep Ar is of the order of 
the collision time of nucleons in a bound nucleus (~ 10"^'^ s). Time intervals are measured by an integer t (number 
of timesteps counted from the beginning of the experiment). 

The real nucleon is simulated by a particle of mass m; this particle is either at rest (symbolized by the zero vector 
Go); or in a state of motion with a single absolute value of the velocity 



(1) 



Accordingly, in our framework the allowed states of motion of a particle are v = iue^, where e,; represents the unit 
vector along the lattice axis i = x,y, z. A CA particle a then exists in one among 7 possible dynamic states Vq, = 0, 
ve^, —ve^^ vey, —vGy, or ve^, —vBz] the momentum of this particle is denoted by (= toVq). As in 19] we do not 
take account of the charge of the nucleons. 

Our CA particles obey an exclusion principle, in the sense that a cell is not allowed to contain more than one 
particle in the same state of motion. Accordingly, the maximum density of our nuclear matter is 7 particles per cell. 



B. Particle Interactions and Dynamics 



To simulate the interactions of a given nucleon a with the rest of the nucleons of our system we choose an 'interaction 
neighbourhood' of cell r containing nucleon a, 7Vint(r). This neighbourhood is the collection of cells made up of the 
central cell r, the tf=& cells which have common faces (0), the e=12 cells which have common edges (e), and the 1^=8 
cells which have common vertices (v) with the 'central' cell r. Nucleon a in cell r interacts with any nucleon in a cell 
r' if and only if r' e iVint(r). 
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The pair interaction energy between a particle in cell r and a particle in cell r', ypair(r, r'), is represented by a step 
potential: Vpair(r,r') = if both cells have a common face; V^air(r, r') = K, and V^air(r, r') = Vy, if they have a 
common edge, or a common vertex respectively. For a pair of particles in the same cell r we adopt a potential of the 
form Vpaii {v,r;p) = Vq + {p — 1)AV; the pair-parameter p takes account of the effect that the interaction energy 
between a pair of particles in the same cell depends on the number of different pairs present in the cell: If there are 
3 particles, and hence 3 pairs, the first pair has an energy Vo, the second pair an energy Vo + AV^, and the third pair 
an energy Vo + 2AV. The total potential in cell r is then the sum of the pair potentials due to all particles in the 
interaction neighbourhood 7Vint(r). 

In order to minimize the number of free parameters of our model we have set = Vo (equal to the interaction 
energy of a single pair in a cell), and V^ = V^ = Vi. We are then left with 3 independent energy parameters specifying 
the pair-interaction potential. The orders of magnitude adopted for the latter are Vo — -3.0 MeV and Vi=-0.3 MeV; 
AV = +1.0 MeV. The precise values are adjusted to obtain optimal agreement with the observations, in particular 
to secure the experimental mean binding energy per nucleon, of -8 MeV in an intermediate-mass nucleus. The value 
of AV has been estimated by requiring that 2 and 3 particles in a cell form a stable bound configuration; a larger 
number of particles per cell leads to an unstable configuration. 

Due to the discrete nature of the allowed CA states of motion, a particle suffers a change of momentum which 
obeys Newton's law of motion in a statistical sense only. Two or more particles in a same cell undergo a scattering 
which satisfies linear momentum conservation. The computational details of the treatment of the transitions among 
the particle states of motion are given in ,1^ . 

C. Fragmentation Clusters and Cluster Configurations 

A major aim of our simulation consists in constructing the distribution law of the fragmentation clusters as actually 
registered by real laboratory detectors, in the form of the number of fragments, N, against size, a (and at fixed time 
t), N = N{a, t). In our experiments the 'size' of a cluster is understood as the number of particles in the cluster. The 
particles are indiscernible, so that any permutation among them which does not alter the occupation of the individual 
cells, does not produce a new cluster. 

A fragmentation cluster of size a of our CA context is eventually identified with a fragment of the laboratory 
nuclear fragmentation process. It simulates a nucleus containing a ~ A nucleons. The experimental counterpart of 
our theoretical distribution, N(a, t), which ignores the charge of the fragments, is then the distribution of the isobars 
as resulting from the laboratory fragmentation process. 

The precise specification of what we understand by a 'cluster' of size a that takes due account of the interactions 
included in our model is given in the Appendix. For the purposes of computing reaction temperatures, we need to 
evaluate the number of distinct configurations ('multiplicities') of the different cluster-geometries compatible with a 
cluster of specified size. This question is also dealt with in the Appendix. 

III. THE SET-UP OF THE NUMERICAL EXPERIMENTS 

All of our experiments are carried out in the CM reference frame whose origin coincides with cell (0, 0, 0). 

A. Initial Configuration 

The initial conditions for a dynamic run are as follows. We simulate two nuclei, referred to as the 'incident nucleus' 
I and the 'target' T, by two identical clusters (in the sense of our definition) located in the half-lattice a; < and 
a; > respectively. The centres of mass of the clusters / and T are required to lie on the x-axis. The mirror-symmetry 
demands that we have rj = {—xm ,0,0) and tt = (+a;Af,0, 0) respectively. The number of particles (and mass, in 
units of the nucleon mass to) of each cluster is At — Aj — A = 150. 

Initially we assign each cluster cell a single particle, so that the initial volumes of the clusters are V — A X^^. The 
shapes of these clusters are chosen to approximate densely packed spheres. This geometry leads to a radius-mass 
relation of the standard type 



/ o \ 1/3 



(2) 
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Overall consistency with the empirical R-A relation then requires A « 1.9 fm. The value we have adopted for our 
model is A = 1.95 fm, which produces the matter density of a real nucleus of mass A = 150. At time t — the two 
clusters are just in contact (Fig. 1). 



B. Initial Motion 



Dynamically, all particles of each initial cluster, / and T, are in ordered motion. In the CM frame of T + /, any 
particle of the target cluster has a velocity — —ve^, while the velocity of any particle in the incident cluster is 
v/ = —vt — +vex. The ordered microscopic motion accounts for the initial macroscopic motion of the two clusters 
along the x-axis. The linear momentum of the global system T + / is zero in this frame. The available impact energy 
is 

E,^p ^ 2 AK , K = ^mv\ (3) 

This energy is transformed into excitation energy of the compound cluster; it leads eventually to the break-up of 
the latter. (As in the traditional statistical models, we ignore here particle-creation processes, in particular pion 
formation). The actual collision occurs at timestep t = 1. Prior to the collision the model describes two clusters T 
and / in uniform motion in the CM frame, with opposite linear momenta = —A mv s-nd p/ = -\-A mv = 
—pT- The velocities of the two clusters are opposite, of absolute value 

|vr| = |v/| = V. (4) 



C. The Cluster Detectors 



In the laboratory experiments the detectors, which identify and count the fragments, and measure their kinetic 
energies, are ideally distributed isotropically around the collision centre. We achieve an acceptable approximation to 
isotropy respecting the lattice symmetry of our CA environment by placing our theoretical counting devices on the 
6 faces (CA lattice planes) of a cube of size 2dD + 1 centred at the origin of the lattice (limiting planes defined by 
X = idjy] y = idjj; z = dzdi)). In conformity with the laboratory experiments, the distance d^ must be chosen 
macroscopically large. The minimum requirement is that do exceeds the distance di„t over which the fragments 
interact (freeze-out radius). For our typical cell size A and for initial clusters of 150 nuclei each moving along the 
X-axis test runs indicate that dint,x ~ dint,y ~ dint.z ~ 20. In our numerical experiments we have chosen djj as large 
as computationally possible {do w 50 for our CA universe of size L=127). 



IV. FRAGMENT IDENTIFICATION AND COUNTING ALGORITHM 



In models of statistical equilibrium, in which the fragments are confined to a fixed finite volume Vb, the cluster 
identification and count can be carried out with a standard algorithm that consists in scanning the whole box available 
to the fragments However, in the real laboratory experiment the detectors are not uniformly distributed over 
a volume, so that the standard algorithm does not duplicate the experimental procedure. The fragment counting 
method we have set up in our simulation is devised to reproduce the principle of the laboratory counting procedure. 

We isolate the 'new' clusters, AAf(a;t;x+), which 'pierce' face x = +dD (our x~^ counter) of the cube at timestep 
t. The total number of clusters of size a, N{a; t; a;+), which have been traced by counter x~^ up to step t, is then the 
sum of all 'new' clusters identified from timestep 1 up to timestep t. We have 

N{a;t;x+) = AN{a;t;x+) + N{a;t-l;x+) 

11 I (5) 

= AN{a;t;x+) + AN{a;t~l;x+) + . . . + AN{a;l;x+) . 

The identification at the other counters x^ , (face x = —djj); y^, (faces y = ±d]j); z^, z^ , (faces z — dzdjj), 
follows the same scheme. The total number of all distinct clusters N{a;t) traced up to time t is the sum over the 
measurements of all 6 detectors. 

The relative distribution of the fragments, Pr{a), probability of a fragment of size a, is given by 

2A 

Pr{a) = £rn with Ar(t) ^ 5^A^(a;0. (6) 

^ ' a—l 
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Finally, if we repeat the same simulation a large number r of times, we have 

where N^^^t) is the total number of fragments of all sizes, N^''\a;t) is the fragment distribution, and Pr^'^\a) is the 
relative fragment distribution in the c*'' experiment. If the number r of experiments is chosen large enough, a stable 
and well defined distribution is expected to emerge. 

To find (Pr(a)) we have developed a straighforward algorithm for the identification of fragmentation clusters of size 
a which enter the plane x = +dD at step t. The algorithm exploits the property that since a free fragment propagates 
with speed v, AN represents the number of fragmentation clusters which intersect the plane x = +dD, but which do 
not intersect the plane x = +d]j + 1. 

To the extent that our CA model incorporates an acceptable approximation to the physics of the fragmentation 
process, the relative fragment distribution as derived from our simulations should duplicate the laboratory distribution 
Pr^^°'^\a) as measured by real detectors. In fact, consider a laboratory nm of total duration t^xp, operating under 
stationary conditions. We then have a constant flux of infalling particles colliding with a flux of target particles. 
Under typical operating conditions these collision processes are binary collisions (each individual collision process 
occurs independently of the other collision processes). Collision c therefore produces in the array of detectors a 
distribution given by Eq. (6). If r{texp) is the total number of collisions which occur over the whole run of the 
experiment, then the detectors register the following relative fragment distribution 



The notations adopted in these expressions are essentially the same as in the theoretical case {N^'^\a;t), iV('^)(t): 
fragment numbers in collision c; Pr^''"'''''^\a): relative fragment distribution of the c*'' collision; Pr^'"^'(a): final 
average relative experimental distribution of the fragments). 



V. DYNAMIC RESULTS 



The evolution of the particles of the two colliding clusters is followed with our CA program over a total time-interval 

twax not exceeding (L — 1) — (dint + do) (about 70 in our experiments). The order of this time interval is fixed by the 
observation that it is the shortest time it takes an individual particle, or a cluster, ejected in the collision, to migrate 
through the available CA lattice, to be reflected on the boundary of the lattice universe, and to be sent back to a 
detector plane. For times t > t^ax a reflected fragment could collide with outflowing fragments; this would violate 
our assumption of non-interaction of the fragments at distances exceeding dint- 
We begin with a discussion of a first physically realistic effect. 

(i) At the collision site we observe a central concentration of nucleons (a compound nucleus) of size Ac{t). This 
concentration progressively loses individual particles as well as fragments. The sequence of frames (a), (b), (c), (d) 
exhibited in Fig. 2 and corresponding to an impact energy per nucleon Eimp/A = 3.75 McV, illustrates this situation 
in greater detail. Frame (a) {t = 20) indicates that besides isolated particles leaving the centre two large fragments 
were symmetrically ejected along the collision line; smaller fragments were blown along the y- and 2;-axes normal to 
the collision line (x-axis). The collision site remains a high density zone, which is clearly visible on the later frames 
(b), (c), (d) {t = 40, 50, and 60). An evaporation of small fragments and individual particles continues from this 
central condensation zone, at a rate which is expected to depend on the impact energy, Ei^p (the only free parameter 
of our series of experiments) . 

The experiment seems to suggest that the collision phenomenon is actually a two-stage process. In the first phase 
of the collision, the nuc;leiis is broken up into several large fragments; the latter are the larger the higher the impact 
energy is; for a high enough impact energy, only 2 fragments emerge from the collision (negligible interaction between 
the two nuclei). The second phase, starting at some timestep to, is a more gentle evaporation from a central fragment. 

To quantify the phenomenon of emission of clusters we fix a small reference volume (a cube of 21'^ cells) which 
encloses the compound nucleus, and which we regard as a rough approximation to the space occupied by the latter. It 
is then reasonable to assume that the rate of particle loss from this reference volume is proportional to some a priori 
unknown power 6 of an excess number of particles over an equilibrium number, Ag, in the reference volume, or 

d/dt Ac = - PiEimp) {Ac - Ao)^ , PiEimp) = 7 Eimp^/\ (7) 
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The functional form of l3{Eimp) follows from a dimensional scaling argument. In our series of experiments in which 
the impact energy alone is regarded as a parameter, all other free coefficients being held fixed, we can express mass 
in terms of the nucleon mass, and volume in the reference volume, so that mass and volume (and hence length) 
become dimensionless; in this system the impact energy has the dimension of an inverse time squared. Hence the 
above relationship. The remaining factor 7 depends on the model constants (i.e., reference volume, mass of nucleons, 
nuclear interaction constants, etc). 

Integration of Eq. (7), from the moment to, at which evaporation starts, to the current time t, yields 

{Ac{t)-Aof-^ = {Acito) - Aof-^ - {1-S) {t-to) Eimp^/^ J, if 5^1, (7.a) 

Ac{t) - Ao = iAcito)-Ao) exp[-^i„//2 7 (i-to)], if S = 1. (7.6) 

Fig. 3 shows the curves {Ac{t) — Ao)'' /{Ac{to) — Aq)'' for values of k (=1- S) ranging from 0.1 to 0.5 (fine lines), 
as well as ln[Ac(t) — ^o]/ln[^c(i) — Ag] (heavy line); time is measured with to taken as origin. The parameters to 
and Ao are read off from the numerical run (to = 11; = 36 in the case of the experiment exhibited in Fig. 3). 
The curves in Fig. 3, which are strongly nonlinear for larger values of fc, tend to become linear in the limit fc — !■ 0, 
as is found from a linear regression analysis (over the 25 timesteps shown). Hence Eq. (7.b) {6 = 1) gives the best 
fit. In other words, the excited compound cluster essentially suffers a standard disintegration. The decay time for the 
specific experiment we have displayed is r^js = 6.58. 

(ii) A second conspicuous feature of the plot of our numerical results is an artefact of the CA lattice symmetry and 
computational procedure. In each lattice direction, zte^, i = x,y, z, we observe a column of particles or clusters that 
propagate all with maximum speed v (Eq. 1) away from the collision site. 

Fig. 2 discloses that the density of particles in each column increases with distance, with a peak density at the two 
ends of each column. The peak density corresponds to the effect of the violent break up occurring in the first phase of 
the collision. The subsequent gentle decrease in time of the content of the residual compound cluster, Ac{t), implies 
that the rate of evaporation (Eq. 7) decreases with time as well. Therefore the density in the columns decreases as 
we approach the central compound cluster; at the latter central cluster, the common centre of the columns, a density 
maximum survives. 

Statistically the particle distribution at any time t preserves the symmetry of the initial configuration. Fig. 2 
demonstrates indeed that the pattern possesses the following symmetry elements: The x- axis is a fourfold axis; the 
planes {x,y) and {x,z), as well as their bisectors are reflection planes; the y and 2;-axes are binary axes, and so are 
the two diagonal axes; (y, z) is a reflection plane. The invariance group of this statistical configuration is D^h (in 
standard notations for the point-groups; cf Landau-Lifshits, Quantum Mechanics) . The distribution of the particles 
along the collision axis x is clearly seen to be different from the distributions along the y or 2;-axes, (cf in particular 
the large fragments on the x axis). The collision axis remains a privileged direction all over the experiment. 

At timestep t = 0, the initial conditions of the collision set-up in our CA lattice environment, Fig. 1, create a 
symmetry breaking from the original octahedral (plus translational) symmetry of the empty (infinite) lattice (the 
'vacuum state'). Oh. to the symmetry of our initial configuration. At all later times, t > 0, our CA results 
demonstrate that the latter symmetry is statistically preserved (cf Fig. 2). The collision process itself, starting at t 
= 1, induces no overall geometric symmetry transition. The occurrence of such a transition would be direct evidence 
for a second order phase transition. If a second order phase transition does occur in the fragmentation process, then 
it must be related to a finer symmetry breaking not immediately manifest in the spatial distribution of the fragments. 

The initial conditions for a head-on collision of two identical nuclei occurring in the continuous physical configuration 
space break the symmetry of the original 'vacuum' (spherical plus translational symmetry) transforming it into a 
cylindrical symmetry, of axis coinciding with the collision axis, Kh Dooh- Provided that our discrete CA model can 
capture the essential physics of the real fragmentation mechanism, the results of the above numerical experiments are 
indicative that the full cylindrical symmetry should be preserved in the real laboratory experiment. This symmetry 
is indeed consistent with the An detector measurements. 

Our next task is an attempt at wiping out the CA artefact (ii) in the computed spatial distribution and to convert 
the latter into a form more directly comparable with the laboratory experiments. 

We reconstruct a Dooft^symmetric number distribution, p'''^^(r; t), from the CA D4fi symmetric number distribution, 
p'-^^{r,t), by the following procedure. We expand the new distribution in Legendre polynomials 

00 

p'-^'=(r;t) = J2 Pnir;t) Pn{cos0), (8) 

n=0 

where the angle is the angle between the position vector r and the collision axis (oriented from the / to the T 
nucleus). Since the CA experiments indicate that the difference in the distribution along the collision axis and the 
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axes normal to the collision axis is relatively small, we truncate the expansion after the dipole terms. In this lowest 
order approximation we then write 

p''^=(r;t) ^ p,{r;t) + r^{r;t) cos0. (8.a) 
The functions pdf", t) and ri{r; t) are then determined from the CA distribution by 

Pcir;t) = B {p^^{r;t))y,,, 

Pcir;t) + r^{r-t) = B (p^^(0, r, 0; t)),, (8.6) 

where the righthand sides are averages over the positions r = \y\ and \z\ on the y and z axes, and r = \x\ on the x 
axis respectively; i? is a normalisation constant which is fixed by requiring that the integration of the distribution 
over the available space (the sum of p'^^'^{r; t) over all CA cells r) is equal to the total number of nucleons (=300 in 
our experiments). 

Panel (b') of Fig. 2 shows a stereoscopic plot of the symmetrized space distribution (8) reconstructed from the CA 
distribution displayed in panel (b). The plot is generated by a straightforward Monte-Carlo procedure distributing 
300 particles in conformity with the statistical law (8. a). We observe that this distribution falls off with distance r 
from the centre (with an exponent < 2), except that the higher density of the distribution front survives. We should 
point out that the method just redistributes the individual particles without conserving the clusters. 

(iii) The cluster-size distribution as generated by our CA dynamics and registered by our counters x"'",etc, obeys a 
power law. This result is illustrated in greater detail in Fig. 4 , which exhibits two instances, one at very low impact 
energy, Eimp/-^ = 0.307 MeV (panel a), and the other at a 13 times higher energy Ei^p/A =3.973 MeV (panel b). 
The plot shows (a non-normalised form of) relation (6) in log-log format, obtained from a total number of runs r = 
16. Over the range a < 10, both curves appear as straight lines, with negative slopes r = 2.56 (low energy) and 2.64 
(high energy). Our experiments indicate that this property holds at least up to Eimp/A « 10 MeV, where r takes 
the value 2.65. We are entitled to conclude that in the range of impact energies Eimp/A < 10 MeV, the (smaller) 
fragments obey a power law with negative slope r « 2.6, which is independent of the energy. As the impact energy 
increases larger clusters are being formed; the tail of the distribution then tends to become longer. 

Critical exponents r of 2.64 and 2.65 have been measured in the laboratory, for the fragment distribution of proton- 
Kr and proton-Xe fragmentation. We should mention also that_in the field of ion cluster fragmentation, fragment 
distributions obeying power laws of exponents of 2.56 and 2.63 ^ have been isolated. 

It appears that the slope values of the fragment distribution we obtain from our dynamic results are in better 
agreement with the laboratory data than the values derived from standard statistical theories, which are close to the 
value given by mean field theory (0|: r = 2.33). This remains also true for ion cluster fragmentation problems, even 
though our model parameters were not adjusted to that particular situation. 

VI. FORMAL TEMPERATURES. COMPARISON WITH LABORATORY EXPERIMENTS 

With the exception of the symmetrised distribution (8), the results described in the previous Section are direct 
results of our A'^-body simulations. They involve no extra approximations besides the assumptions inherent in any CA 
modellisation (space and time discretisation) and the schematisation of the interaction potential among particles. In 
particular, the dynamic calculations dispense with the hypothesis of formation of a compound cluster in the collision 
process. Our numerical experiments indicate that a configuration showing a higher central concentration, which we 
refer to as the compound cluster, always emerges in the simulation (Fig. 2), provided only that the collision energy is 
low enough as compared to the binding energy of the nucleons in the nuclei. The compound cluster (i.e. high central 
concentration) survives over a period of time which can be estimated from Eqs. (7). 

Within this compound cluster, and, in principle, also within our CA universe U {127^ cells; toroidal topology), a 
variety of statistical equilibria e, e', e", . . . are conceivable among the different components. 

(a) If we identify the components as the individual particles (nucleons in our specific nuclear fragmentation simula- 
tion; etc) of a configuration then we may have a statistical equilibrium in Z, ce^^'^^K due to energy exchange (E) 
among the particles {n{Z)) of configuration Z. The configuration Z may be the compound cluster C, or our entire 
CA universe U. 

(b) If we identify the components as the clusters of particles c (which in turn, when evaporating from the compound 
cluster become the observable fragments) , then we may have (i) a statistical equilibrium due to energy exchange among 
the clusters of a configuration Z, c{Z), ce'^^^'^ ■ But we may also have (ii) an equilibrium due to exchange of energy as 
well as particles (P) among these clusters, e^j^p'^^'^-'. Under the latter alternative, the observed distribution of frag- 
ments against size will be the thermal reaction-equilibrium distribution of the clusters. For the sake of completeness 
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we mention further (iii) that within each individual cluster c/, the component particles may exchange energy among 
each other; this process may then lead to a special type of equilibrium, eE"'^'^'\ inside each cluster c/. 

We define here a temperature parameter T(e), which is associated with a specific thermodynamic equilibrium 
process e, as follows : The parameter T(e) is the thermodynamic temperature that reproduces the macroscopic 
properties of the specific equilibrium e assuming it is realized. 11 uj = F{e;T) is the thermodynamic relation that 
expresses the observable w as a function of the thermodynamic temperature T under equilibrium conditions e (other 
thermodynamic variables being held fixed), then from the measurement of uj and the thermodynamic relation we 
set r(e) = T. (Note that our way of introducing temperatures differs from the more formal procedure adopted in 
Chernomeretz et al. |2^). We are free to use this thermodynamic relation in a formal way to determine a parameter 
T(e) characterising a given macroscopic configuration, whether or not the specific equilibrium e holds. We begin with 
listing the thermodynamic relations of relevance for our purposes. 

In principle, all of the above-listed equilibria, e, e', e", ... may arise in our numerical experiments, or in the 
laboratory experiments (in the sense that there is for instance no membrane surrounding a fragment that would 
prevent the exchange of particles, etc). The question is rather: Given an equilibrium process e, does the collision 
system we investigate survive over a time-span that is long enough for the equilibrium e to establish itself? 

If the collision system reaches a full thermodynamic equilibrium, then we must have 

T(e) = r(e') = . . . = T , (9) 

where the formal temperatures T(e), . . . obtained as indicated; all of these formal parameters are then equal; they 
define the single thermodynamic temperature T. Conversely, if relation (9) is violated, then the relevant statistical 
equilibria are not realized. We have indeed good reasons to believe that some of the possible equilibria listed above 
will never materialize in our system (cf below); or, alternatively, some equilibria cannot establish themselves over 
certain ranges of the impact energy. Our numerical results confirm this point. 

Even if a state of full thermal equilibrium is not attained for the system we are investigating, the formal tempera- 
tures, T{e), T'(e'), . . ., may remain perfectly useful parameters for the purposes of comparison of the numerical results 
with actual laboratory experiments in the following sense. Suppose a given observable fl, measured in the laboratory, 
is plotted against the temperature parameter, T(e), = Q{T(e)), this temperature being measured according to a 
well-defined protocol (cf above). A typical instance is provided by the (formal) caloric curve, T{eE,p'^'''^^) vs Eimp, 
(equilibrium process: energy and particle exchange among clusters in the compound cluster). The necessary condition 
for a given theoretical model, such as our present CA model, to be an adequate model for the fragmentation process, 
is then that this model can duplicate the experimental f2 — Q{T{e)) relation. The agreement must hold provided 
that the formal parameter T(e) be obtained in conformity with the experimental protocol. It should be kept in mind 
that there is no guarantee that an easily measurable temperature parameter of our CA gas (such as r(eE"('^)) can 
be substituted to the actual experimental parameter (temperature associated with ratio of yields T{eE.p'^^'^'^) )■ 

We have measured the following collection of temperature parameters for our CA system (adapted to nuclear 
fragmentation) . 

A. Nucleon-gas temperatures: Global nucleon— gas temperature, T„^u){t), and nucleon— gas temperature in 

compound nucleus, T„^c){'t) 

The temperature parameters Tj^(jj'^[t) and Tn(c){t) are the formal temperatures of a gas of nucleons in thermal 
equilibrium under exchange of energy in the CA universe U and in the compound nucleus C respectively. With the 
above notations Tn(u) (t) = 2^(e£;"^'^-' ; t) , T'n(c) (t) = T{e b"^'^-' ; t) ; the extra argument t indicates that these parameters 
depend on time. 

The total number of nucleons of our system, the CA universe U, is conserved, Au{t) = A (=300). The number 
of nucleons in the compound nucleus C at timestep t, Ac{t), is variable, and so are the number of nucleons in 
motion in the universe, ^m(c/)(^)i and the number of nucleons in motion in the compound nucleus, Arn(c)(^)- The 
temperature parameters, which are measures of the average kinetic energy per nucleon of the gas of nucleons in U 
and C respectively, are then given in terms of these numbers of nucleons by the following expressions 

I kT^iu)it) = 2_^J2M)r = ^^^K < K, (10) 

a— 1 

2 "(ojv / 2 Ac{t) ^" Acit) ^ ' 
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In these expressions k is the Boltzmann constant, Va{t) is the velocity of the nucleon labelled a and taken at 
timestep t, and K is the kinetic energy of a moving CA particle (Eq. 3). Note that all formal temperature parameters 
we shall introduce are time-dependent. 

The second temperature parameter, r^^p) (Eq. 11), measures the physically meaningful average kinetic energy of 
a gas of particles in the compound cluster, which can in principle reach a thermal equilibrium in our CA model, or 
in the laboratory; hence it can represent a genuine temperature of a gas of nucleons. This is the case if the energy 
exchanges among the nucleons have time to establish themselves, i.e., if the lifetime of the central condensation (cf 
Eqs. 7) exceeds the average collision time among nucleons in C (a few timesteps). 

Figure 5 shows a caloric curve Tn(c) against the impact energy per nucleon Eimp/A, with Tn(c) = Tn{c}(tref) 
computed from our CA runs at a reference timestep tref (= 17 for reasons to be discussed below). In the low energy 
range, up to E^/A « 7-8 MeV/nucleon, the temperature parameter rises slowly with impact energy, approximately 
as 

'^n(C) = 0.15 Eijnp/A, Eijnp < E^- (11. a) 

In the high energy range the rise is steeper, and we approach 

2 

Tn(C) = 2 {Eimp ~ Ec)/A, Eimp ^ Ec- (11-^) 

The slope 2/3 in the latter relation is indicative that the random kinetic energy of the nucleons is asymptotically 
equal to the impact energy minus E^- Accordingly this critical energy appears as a binding energy of the nucleons in 
the compound nucleus. The numerical value of Ec/A, which is of the order of the average binding energy per nucleon 
in a stable nucleus, is indeed consistent with this interpretation. 

As can be seen in Fig. 5, the slope of the asymptotic expression (ll.b) is in line with the high energy branch of 
the experimental caloric curve of the ^^^Au-^^^Au fragmentation process [lj|. The NuPECC caloric curve exhibits a 
plateau, at a temperature level of 4.5-5.0 MeV, which is not reproduced by the CA nucleon-gas temperature; however, 
the change of slope in the theoretical curve coincides with the transition from a plateau to a rising behaviour. We 
should stress that the experimental points of Trautmann's analysis, also plotted in Fig. 5, indicate no proper plateau; 
they rather follow a rising curve similar to the CA Tn(c) curve. Quantitatively the experimental data (NuPECC or 
Trautmann) are shifted by AT w 3.5 Mev above the CA curve. 

The first parameter, T^^u) must be regarded as an artificial temperature. It measures an average kinetic energy of 
nucleons composing a mixture of two gases, which in the context of our CA model, or the real laboratory experiments, 
never interact. Namely we have on the one hand a gas of nucleons essentially trapped inside the compound nucleus 
C, and on the other hand an expanding gas outside the compound nucleus, U — C . In the latter gas the collisions 
are negligible; no energy exchange can take place among the outer nucleons. All of these latter nucleons conserve the 
momentum and kinetic energy (maximum energy K) they have acquired at the moment they evaporate from C . To 
simulate realistic laboratory nuclear fragmentation, our CA experiments must consist in relatively short runs, of a 
total number of timesteps tmax essentially chosen as follows (cf Section 4) . Once a fragment has left the collision site 
it suffers no further interactions. In the finite CA universe of our model this requires that an ejected fragment be not 
allowed to be reflected on the boundaries of the cubic CA universe, and sent back to the reaction site. Very roughly 
we then choose tmax < L/2. (Interactions in [/ — C would require that the number of timesteps obey tmax ^ L). 

Even though T,i(ij) (0 has a formal meaning only in the nuclear fragmentation problem, it can be used to find the 
physically significant temperature parameter of the compound nucleus, r„(p')(t). In fact, in the numerical experiments 
we can easily compute Tnijj){t)] we can also easily count the outer nucleons, Ajj^c{t)- Taking then account of the 
relations 

A = Ac{t) + Au-c{t), and Am(u){t) = A«(c)(t) + Au^c{t), 
we obtain from Eqs. (10,11) 



= r„,„(i) - (H f - T„,„,«)) -^^^y (12) 

This relation demonstrates that the temperature parameter in the compound nucleus at timestep t, r„((7)(i), is always 
less than the formal temperature parameter of the CA universe, T'„([/)(t), taken at the same timestep t. Fig. 5 exhibits 
also the formal caloric curve Tn(u){t) vs Eimp/A, at timestep t — tref (= 17), which is seen to obey the inequality 
Tn(u){t) > Tn(^c){t)- In the lower energy range the two curves Tn{u){t) and Tn(c){t) are nearly superposed, indicating 
that at the reference time few nucleons have escaped the compound nucleus. Even at the highest impact energies we 
have investigated, the difference in the two formal temperatures does not exceed 1.5 MeV. 
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B. Cluster-gas temperatures Tc{ij){'t) and Tc(e)(i) 

At timestep t consider the specific /"^ fragmentation cluster, c/, / = 1,2,. . .,Fz{t), a representative of the cluster 
equivalence class C{af). Fz denotes the total number of clusters in the configuration ([/, or C). Denote by (vj(t)) 
the velocity of this cluster (average of the velocities of the component nucleons) at step t. Then the cluster-gas 
temperature is the temperature parameter associated with the equilibrium brought about by the exchange of energy 
among the clusters in a given configuration Z. We have Tc(z)it) = F{eE'^^^'^',t), Z—{U,C). Hence 

3 1 

As in the case of the nucleon-gas, the cluster-gas temperature parameter (t) is again a formal magnitude, 
since the cluster-gas outside the compound nucleus does not interact within our runs of t^ax timesteps. On the 
other hand the clusters in the compound nucleus C may have time to thermalise, so that the temperature parameter 
Fc(c)ii) does represent a physically meaningful temperature. 

On Fig. 5 we have superposed the formal theoretical caloric curves Tc([/)(t) and rc((7)(i) vs impact energy Eimp/A 
(at timestep ire/)- We note that the inequality between the two temperature parameters of the nucleon gas in U and 
C, read off from Eq. (12), is preserved in the case of the cluster gas in U and C : Tc^u^{t) > T^(^c){i) (outer fragments 
move with maximum speed). The cluster-gas curves have much steeper slopes than the nucleon-gas curves, implying 
the further inequality Tt,(c)(t) > T„(c)(t). 

The higher cluster-gas temperature reflects the following property. In a given configuration Z the number of 
clusters is smaller than the number of nucleons while the total energy to be shared among clusters, or among nucleons 
is the same. Therefore the average energy per cluster is larger than the average energy per nucleon. At high impact 
energies, Ei^p/A > average binding energy per nucleon in a stable nucleus, one might expect intuitively that the 
physically meaningful caloric curve for a fragmentation process should be the nucleon-gas curve Tn(c)- The high 
energy available, when shared among the nucleons, would indeed produce average energies per nucleon which exceed 
the binding energies of any fragment. However, our CA experiments demonstrate that in the case of high impact 
energy the merged T and / nuclei break immediately up into essentially 7 large fragments. Two compact fragments 
carrying a sizeable fraction of the mass of the / and T nuclei, continue to propagate along the collision axis. Four 
smaller fragments are ejected along the y- and z-axes respectively. The residual nucleons form a concentration of 
matter at the centre. Only the latter can play the part of a compound nucleus in which actual energy sharing may 
occur. Qualitatively this initial break-up survives at lower energies, except that the 6 ejected fragments become 
progressively smaller while the central residual becomes larger as the impact energy is decreased (cf the sequence 
shown in Fig. 2 for an impact energy 3.75 MeV per nucleon, where the 7 fragments are clearly distinguishable on all 
panels). 

The observed scenario is indicative that at high impact energies it is the cluster-gas temperature that supplies 
the physically meaningful characterization of the laboratory caloric curve. This temperature takes properly care of 
the contribution of the large fragments in the energy balance. Any laboratory measurement technique, (and any 
theoretical procedure) of a temperature assignment ignoring the largest clusters must fail to provide a physically 
relevant temperature estimate. 

This point is of importance in connection with the reaction temperatures. The latter, to the extent that they 
typically refer to fragments of low size (cf below), are not physically representative in the high energy range. 

C. Reaction temperatures, Tcci...id' ...(^) 

The favoured laboratory method for assigning an experimental temperature to a fragmenting nuclear system con- 
sists in measuring ratios of yields of different fragments, assuming a statistical equilibrium of general type eB^p"^*-*^^ 
among the fragments [231 . Examples are the ratios ('^He/'^He) and (^Li/^Li) for the earlier analysis of the ALADIN 
experiments Il3l . [2^ a nd , and other ratios of populations of isotopes in the more recent analysis of the same 
experiment llll; (cf also '24^). In "2^ the ratios (3He/4He)/(d/t) and (3He/4He)/(^Li/7Li) are considered for the 
multifragmentation resulting from an Au target bombarded by C ions. Several other instances of measured ratios are 
listed in 0. In all cases populations of light nuclei alone have been investigated, from which a specific temperature 
parameter is derived via the standard relation for chemical equilibria (cf ; the specific nuclear context is dealt with 
in H). 

In the framework of our CA formulation which ignores electric charge, the chemical equilibria among different 
fragments are not directly comparable with the real nuclear equilibria investigated, which involve isotopes. Under 
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conditions of true thermal equilibria, the ratio among any group of isotopes, or of isobars, is governed by the same 
thermodynamic temperature. But if the system is not in a genuine state of equilibrium under nucleon exchange, 
then the temperature parameter is specific for the precise reaction process. Different reactions, and hence different 
measured ratios of fragments in the CA experiments, lead to different temperature parameters. The same conclusion 
holds for the laboratory ratio measurements. 

In fact, as transpires from |l4j . different experimental ratios, and different analyses of these ratios, have led to 
different shapes of the caloric curve in the case of the ^^^Au-^^^Au fragmentation process. In Trautmann Il3| a rising 
pattern for the He-Li ratio is identified, while a nearly constant temperature is found for other ratios; moreover, the 
earlier NuPECC caloric curve based on other ratios, exhibits a temperature-plateau (not present in Trautmann), 
approximately over the range 3 to 10 MeV per nucleon. In Fig. 5 we have superposed the various available experimental 
points defining caloric curves for the symmetric -'^^^Au-^^'^Au fragmentation. This laboratory process comes close to 
our CA simulation, even though the total number of nucleons involved in the laboratory is 4/3 times the number of 
nucleons of our simulation. 

To formulate the relevant statistical expressions for reaction equilibria in the CA context, consider the equilibrium 
among the cluster classes C(a/), C(a/'), . . ., described by the stoichiometric scheme 

lyf C{af) + Vf, C{af,) + ... ^ Vf,, C{afn) + (14) 

{uf, Vfi,. . . integers consistent with conservation of the number of nucleons in the reaction process : Vfaf + VfCLfi + 
• ■ ■ = ^/"fl/" + ■ ■ •) 

Denote by N{af,t) the total number of clusters of class C{af) present in our system at timestep t. If a statistical 
equilibrium holds, then the standard statistical procedure allows us to write 

V(2TrmkT(t)f/^ af^/^ ( , , ^^ \ 

^(«/'*) = ^ j;r-^ ^ ^ 1 ^ eM-E^nt.J/kT{t)) 1 X exp(Aiy/ a/). (15) 

In this expression V denotes the reaction volume; T{t) is the equilibrium temperature at timestep t. The summation 
extends over the J internal energy states Eintj of the J distinct geometric configurations of the same cluster class 
C{af) (as defined in the Appendix). Physically these energies simulate different excitation states of a nuclear fragment 
of mass number aj. The factor gj is the statistical weight of the energy state Eint,j- 

To compute the parameters referring to the internal states we construct the different geometrically different cluster 
classes denoted G(a/,c/;7) in the Appendix. We then evaluate the corresponding energies, E{af , Cf;"f), and the 
related multiplicities g{af , Cf;"f). The full details are given in the Appendix. The factor A is the Lagrange multiplier 
that takes account of conservation of nucleons. The factor is the 'volume' of an elementary phase-space cell. In the 
strict classical context of statistical mechanics the elementary cell is not defined; a quasi-classical argument identifies 
h with Planck's constant j^^. In the CA context we have a natural phase space cell, inherited from the discretized 
space and discretized velocity, or momentum space; the 'volume' of this cell is (AA x m x v)'^ which is to be substituted 
to h^. 

We re-write Eq. (15) in the form 

N{af,t) = Q{T{t),V)xa{af)x Z,nt{af,T{t))xeMXvf af). (15.a) 

The first factor of the righthand-side, Q{T{t), V) is essentially the translational partition function of the clusters of 
class C{af) (normalised to a/ = l; the actual mass contribution, af^^'^, is included in the second factor a{af)). The 
translation contribution is evaluated in the standard context of classical mechanics, for reasons of algebraic simplicity 
(the discrete kinetic energy states of the CA lead to a more complicated expression, which should be equivalent to 
the classical expression in the large a/ limit). This factor is the same for all species of clusters. The second factor, 
cr{af) = aj^l"^ l(aj — 1)!, is a combination of the mass effect in the kinetic energy contribution, and the effect of 
the indiscernability of the nucleons (invariance under permutation of all nucleons of the cluster). The third factor is 
essentially the internal partition function of the cluster class C(a/), which we evaluate in the specific CA framework. 

For the purposes of estimating reaction-equilibrium temperatures from our numerical experiments we restrict 
ourselves to clusters of smallest sizes, a/ = 1, 2 and 3. As transpires from the Appendix, the precise enumeration 
of the cluster geometries becomes already rather involved for a/ = 3. To deal with higher sizes, we believe that 
asymptotic approximations to the cluster configurations should be constructed. This has not been done in the present 
work. 

o-S = i 

For a cluster made of a single nucleon we have (7(1) — 1, and the internal partition function reduces to 



Z,^t(\,T(t)) = 1. 



(16) 



13 



For a cluster made of 2 nucleons (j(2)=-\/2; the internal partition function involves the contributions of the geometric 
configurations listed under Eqs. A3-A5, and A6. 

Zi„t(2,r(i)) = % exp[-y^/fcr(t)] + \ cxp[-K./fcT(t)] 

^ y ^ (17) 

+ - exp[-K//cT(t)] + 1 exp[-K/A;T(i)] . 

«/ =5 

A cluster of 3 nucleons has (t(3)=\/3/2. The geometric configurations which contribute to the partition function 
are listed in the Appendix under Eqs. (A7)-(A9) (first line of Eq. (18)), (A10)-(A15) (second and third lines), 
(A16)-(A19) (fourth and fifth lines), (A20)-(A22) (sixth line), and (A23) (last line): 

^mt(3,T(t)) =1 exph2F0/fcT(i)] + | exp[-2ye/A:T(t)] + ^ exph2K/fcT(t)] + 

((^2/2 -<^) exp[-(2y0 + V;)/fcr(i)] + (eV6-e) exp[-2ye/fcr(*)] + 
eV6 exp[-2V;/A:T(t)] + eV6 cxp[-(2Ve + K)/fcr(t)] + e exp[-2K/A;T(t)] + 

2^. exp[-2K/A:r(i)] + ^ cxp[-(T/^ + K.) Ar(t)] + 
^ exp[-(y^ + + K)/fcT(i)] + ^ exp[-(y^ + K)/A;T(i)] + 

^exp[-(ye + K,)AT(t)] exp[-(2y^ + K)/fcT(t)] + 

e exp[-(214 + V;)/fcT(i)] + v exp[-(2K + Vo)lkT{t)\ + 1 exp[-3(K, + AV)/A;T(t)] . (18) 

The following ratios of cluster numbers are independent of the Lagrange parameter A: 

^(3,t)iV(l,t) ^ V3 ^int(3,r(t)) 

iV(2,f)2 4 Zi„*(2,r(t))2' ^ ^ 

iV(3,i) _ x/3 ^ 1 Z,„t(3,r(0) ^ ^2o.a) 



7V(2,i)Af(l,t) 2^/2 e(T(t),l/) Zi„t(2, T(i)) ' 



^^^^ = X ^.^i X ^i„t(2,T(t)). (20.6) 



The first ratio is independent of the translational term Q (and hence independent of the reaction volume V). The 
numbers of clusters N{af,t), a/=l,2,3 are directly supplied by our CA runs, at every timestep t, so that the lefthand- 
sides of Eqs. (19, 20. a, 20. b) are known. These equations are then solved with respect to the parameter T{t) (written 
Ti3/22{t), T^/i2{t), T2fii{t) respectively, under the above alternatives). These temperature parameters are the reaction 
temperatures of the equilibrium processes 

C(l) + C(3) ^ 2 C{2) , C(3) ^ C(l) + C(2) , C(2) ^ 2 C(l) , 

respectively. 

In order to follow as closely as possible the actual laboratory experiment, consider the cmiiulated number of clusters 
of size ttf, M{af, t), which have passed the collection of detectors up to timestep t. The relevant timestep is identified 
with the length of our run, tmax ~ 70 steps. With our detectors situated at a finite distance do (~ 50) from the 
reaction centre, the first clusters arriving at the detectors were emitted from the collision site do — Rc timesteps 
prior to the arrival time {Rc, radius of the initial merged T + 1 configuration, « 10 in our set-up). Accordingly, in a 
run of 70 steps the first fragments have time to be reflected on the boundary of our CA universe. The actual choice 
of tmax secures that the clusters cannot reach the detectors after reflection; thereby the clusters arriving close to the 
detectors cannot undergo any interactions with other clusters. 

On the other hand, the last fragments registered by the detectors, at step tmax have left the central nucleus roughly 
at step tmax — {do — Rc) (~ 30). Accordingly, the cumulated number readings of the detectors, terminating at 
timestep tmax (=70), M{af,tmax), cover the first 30 timesteps of the fragmentation mechanism. Alternatively, the 
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cumulated number may be interpreted as the average number of clusters of size a / actually present in our system at 
the 'average' time 1/2 [tmax — {do — Rc)] = tref, of the order of 15: 

N{af,Uef) = M{af,t„,ax)- (21) 

In our numerical simulations we have set this reference time tref equal to 17. 

The three formal caloric curves Ti^/22{t)~Eimp / A, T^/ i2{t)~Eimp / A, and T2/ii{t) ~Eimp/A, for t = tref, are plotted 
in Fig. 5. All three curves are seen to be essentially independent of the energy over the range E/A <25 MeV 
investigated: 

T13/22 « T3/21 « Tyn « 1 - 2 MeV. 

Qualitatively this behaviour is in line with the different formal temperatures derived from the laboratory Au-Au 
fragmentation as analysed in Trautmann [3 | (with the exception of the He-Li isotope temperature) . Quantitatively 
the constant temperature level as found in the laboratory experiments lies at 5 Mev. The experimental temperature 
is thus shifted by AT « 3.5-4 Mev with respect to the CA temperature: 

j.Lab ^ J.CA _^ ^22) 

In the case of the nucleon-gas temperature we have noticed a shift of a similar order between the theoretical CA 
temperature and the laboratory estimate. 

The temperature defect AT between the CA and experimental temperature is thought to be due to a feature inherent 
to our CA treatment. The construction of a stable nucleus in the CA framework relies on a discretized version of 
classical mechanics, in which the component nucleons possess no kinetic energy at all (in the reference system attached 
with the centre of mass of the nucleus). Since classically temperatures are related to microscopic kinetic energies, an 
underestimated kinetic energy leads to underestimating the temperature as well. The reaction temperatures encoded 
in relations (19), (20. a), (20. b) refer to small-size fragments only. According to our remarks on the cluster-gas 
temperature, these temperatures are not thought to be representative for the real collection of fragments at high 
energies {Eimp/A > 8 MeV). However, in the low energy domain the available energy can concentrate on small-size 
fragments, which then form and dissolve easily; therefore the latter fragment distribution does reflect a physically 
meaningful temperature at lower energies. The full physically significant temperature run with energy (caloric curve), 
from to about 25 MeV, is suggested therefore to be made of the cluster temperature at the high energy end, 

and the reaction temperatures (T2/11, etc) at the low energy end, with a continuous transition from one curve to the 
other around the critical energy Ec/A « 7-8 McV/nucleon, of the order of the average binding energy. 

VII. CONCLUSION 

The primary aim of the proposed CA simulation has been to devise a framework capable of replicating the actual 
laboratory procedure of monitoring a fragmentation process generated in cluster collisions. Previous theoretical work, 
investigating thermodynamic properties of the collection of fragments (e.g., a caloric curve), was implicitly based on 
an assumption of a thermodynamic equilibrium. In the present model no equilibrium hypothesis is needed. Statistical 
relations (Eqs 10, 11; 13; 19, 20. a, b) are used formally, for the purpose of making comparisons with the laboratory 
experiments. 

Our numerical experiments demonstrate that the distribution of the CA clusters against cluster size, as registered 
by a collection of detectors surrounding the collision site, obeys a power law of exponent r close to 2.6. This model 
result is in excellent agreement with the laboratory results of nuclear as well as other multifragmentation processes. 
It is well known that the very existence of a power law can be derived in the context of various statistical models 
(cf the review 0). The statistical assumptions of these theoretical approaches, and the related counting procedure, 
do not respect, however, the real laboratory protocol. It should then not come as a surprise that the theoretically 
derived slope of the (equilibrium) distribution is found to be significantly different from the experimentally measured 
slope of the (far from equilibrium) distribution (2.3 against thexperimental value 2.6). 

Secondly, the CA experiments are indicative that the notion of a compound cluster, which would support a mean- 
ingful thermodynamic treatment, is of limited value in the fragmentation problem. We observe that typically during 
the earliest phases of the collision process the combined T + / is fractured into a few large fragments. The latter tend 
to acquire all the available mass when the impact energy becomes large enough; in the latter limit, Eimp 3> binding 
energy, two fragments survive. The target T then becomes transparent to the incident cluster I. Significant energy 
sharing is found to occur in the range of low impact energies. 

Thirdly, our CA model experiments demonstrate that the original spatial symmetry of the collision set-up is 
statistically preserved during the whole fragmentation mechanism (Fig. 2). 
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Finally, a major goal of our Paper was to construct formal caloric curves, T^^ - E^^^IA, based on a variety of 
formal temperature parameters directly derived from the CA experiments, and defined and discussed in this Paper 
for clusters of particles of arbitrary nature. In the nuclear fragmentation set-up, comparison of our CA caloric curves 
(temperature parameter chosen at a selected timestep tref), with the laboratory caloric curves, T^^^-E.^np, derived 
from the Au-Au collision experiments in the case of different methods of measurement of a temperature T^*^^, calls 
for several comments. 

With none of the formal temperatures introduced we can reproduce the qualitative shape of the full NuPECC caloric 
curve 14]. The theoretical treatment does not reveal a transition from a rising behaviour to a plateau, and again from 
a plateau to a rising branch. Full CA caloric curves either reduce to a plateau (case of the reaction temperatures, 
Eqs. (19), (20. a), (23. b)); or else the curves are rising everywhere (nucleon gas, Eqs. (10), (11); cluster gas, Eq. 
(13)). The formal CA temperatures are qualitatively more in line with the Trautmann analysis of the experimental 
Au-Au. Experimental reaction temperatures are found to be independent of the excitation energy per nucleon; this 
behaviour is duplicated for our 3 CA reaction temperatures. The laboratory He-Li temperature parameter exhibits 
a rising behaviour, reminiscent of the rising pattern of the temperatures of the nucleon gas and the cluster gas. The 
slope of the latter experimental curve, of approximately 0.27 (in the energy range to 15 MeV), compares favourably 
with the average slope of the cluster-gas, of 0.3, in the energy range < 7.5 MeV. Quantitatively the CA temperatures 
are typically too low as compared to the laboratory measurements of the temperature parameters, by an amount AT 
of nearly 4 MeV (Eq. 22). We have traced this effect to the treatment of the dynamics of the nucleons in a context of 
classical mechanics; in this framework the residual quantum mechanical zero-point kinetic energy is disregarded. We 
believe that this effect may account for a defect in the theoretical temperatures. An extension of the model taking 
account of quantum-mechanical effects will eventually be needed to handle the energy problem adequately (cf (3Q | for 
an attempt at implementing quantum mechanics in the CA framework). 

In the specific nuclear multifragmentation case we have analysed, besides the absence of quantum corrections, our 
model ignores any charge-related effects, so that the detailed theoretical results cannot be compared directly with 
the detailed experimental measurements. All laboratory cluster identifications rely on counts of isotopes (rather than 
isobars, as done in our approach). We hope to be able to include electrostatic effects at a later stage. We also plan to 
extend the model to handle asymmetric collisions {Ai ^ At) and collisions involving a non-zero impact parameter. 

The CA model as developed in this Paper remains a manifestly highly schematic representation of an actual physical 
fragmentation problem of any nature, and as such it can only replicate, and hence also isolate, properties which are 
largely insensitive to the microscopic details of the physics. Our numerical experiments suggest that the power law of 
the fragments, and on a quantitative level, the exponent of the latter, belong into this category of invariant properties. 
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IX. APPENDIX : ENUMERATION OF CLUSTER CONFIGURATIONS 
A. Clusters. Cluster equivalence classes. Cluster geometries 

We introduce a definition of a 'cluster' of size a that rests on the notion of 'interaction neighbourhood' A'int(r) [l^ . 
If r labels an arbitrary cell, then any cell r' distinct from r and contained in the interaction neighbourhood of cell r, 
iVint(r), is termed 'adjacent' to cell r. The notion of adjacency is naturally extended to an arbitrary set of cells, S: A 
cell r' is adjacent to the set of cells S if (i) r' ^ S; and (ii) there is a cell r G S such that r' G Ni^t{r). The collection 
of all cells adjacent to S is the outer boundary of S, dS. 

Define a 'walk' in the CA, of head Vo and tail r/, w{ro,rf), as an ordered collection of cells (ro, ri, . . . , r/_i, r/), 
such that for any pair of successive cells (r^, rj+i) we have r^+i) g iVint(rj) (cf [24] for the graph-theoretical details). 
Define further a 'connected set of cells', Q, in the lattice space of the CA as a set of cells {ra,rf,, . . .}, such that for 
any pair of cells of the set, r„, r„i G Q, there exists a walk w(r„, r„i) € Q (all cells of the walk lie in the connected 
set). 

We understand by a 'cluster rooted at particle a', F[a], the connected set of cells such that (i) particle a is located 
in one cell of the set; (ii) each cell of the set is non-empty (it contains at least one particle); and (iii) the outer 
boundary of the connected set, 9-F[q;], is empty. 

A cluster rooted at particle /3, F[f3], is identical with the cluster rooted at particle a, F[a], F[/3] = F[a], if particle 
/3 occupies a cell of the connected set of cells F[a]. Accordingly we can talk about 'cluster F\ without reference to 
a root particle. A 'cluster', or 'fragmentation cluster', F, is a connected set of non-empty cells, which has an empty 
boundary. 

We distinguish the fragmentation clusters occurring in a CA experiment by an argument /, writing F(f) to refer 
to the 'fragmentation cluster' in a given CA context (at a given timestep t). Any collection of clusters F{i), F{j), 
. . ., containing the same number of particles, regarded as indiscernible, ai = aj = . . . = af, belongs into the same 
'cluster equivalence class' C{af). It is the cluster equivalence class C{af) which we regard as the CA counterpart of 
the 'nuclear fragment' of mass number A — aj in the nuclear fragmentation laboratory experiment. 

A pair of clusters F{i), F{j), of same cluster equivalence class, in which the cells are joined according to a same 
geometric rule, and such that each cluster of the collection has same binding energy, will be said to have same 'cluster 
geometry' G. The collection of all clusters (of same cluster equivalence class) which have same cluster geometry define 
the 'cluster-geometric equivalence class' G. 

Two clusters F{i), F{j) of same cluster geometry G are said to be 'geometrically equal' if after translation along 
the lattice axes, and permutation of the particles among the cells, they can be superposed exactly. Otherwise the 
two configurations are geometrically unequal or distinct. The number g of geometrically distinct configurations of a 
cluster-geometric equivalence class G is the 'multiplicity' of the cluster geometry. 

The reason for assigning a special status to clusters superposable under translations is that the statistical treatment 
of Section 5 deals separately with the translational motions of the fragments, the effect on the reaction equilibrium 
of the energy attached to these motions being described by the function 8(r, V) (Eq. 15. a) (a translational partition 
function). The evaluation of the internal partition function Zint, which concerns us here, involves those configurations 
which we have referred to as 'geometrically distinct'. 

B. Construction of cluster geometries 

Essentially, the notion of cluster geometry enables us to group together clusters which are geometrically distinct, but 
which become superposable after application of certain geometric transformation groups (the discrete lattice symme- 
tries). It is manifest that clusters of same cluster-geometry have generically same binding energies. Conversely, if the 
interaction energies are generic, then different cluster geometries realize different binding energies (accidental degen- 
eracies may occur as a result of a non-generic interaction potential). For the purposes of the statistical mechanics of 
Section 5 'energy equivalence classes' rather than cluster-geometric equivalence classes have to be isolated. The above 
comment indicates that the two questions are essentially equivalent (or closely related in the case of degeneracies). 

We have the following natural construction procedure for a cluster geometry, G(a, c; 7), which we specify (i) by the 
number of particles a in the cluster; (ii) the number of cells c of the cluster; (iii) the precise mode of assembling the 
cells, and the precise distribution of the a particles among the c cells, which we symbolize by the descriptive parameter 
7: 

(1) Assemble the c different cells according to the selected construction rule 7 to form a connected set. 

(2) Distribute the a (> c) particles among the c cells; this procedure generates one representative of the cluster- 
geometry G(a, c; 7). 
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(3) To find the multiplicity g{a, c; 7) of the cluster geometry, generate all geometrically distinct configurations 
representative of the same cluster geometry; this is done by applying the symmetry operations of the lattice (barring 
translations along the axes, as well as those geometric symmetries which are equivalent to particle permutations). The 
binding energy of a representative of this cluster geometry, naturally denoted by E{a,c;^), is obtained by applying 
the rules of Section 2. 

The totality of cluster-geometries corresponding to a cluster C{a) is finally generated by repeating steps (1), (2) 
and (3) for all allowed choices of cells, c = a, a — 1, . . ., 1, and all possible geometrical assemblages of the cells into 
clusters. 

For a/ = 1 there is only one cluster equivalence class; the binding energy is zero. For any af > 1 the cluster 
equivalence class of a specied cluster C{af) contains several cluster geometries G{af,Cf;j), which are energetically 
distinct. 

In order to carry out the construction and characterization of the cluster geometries it is helpful to resort to an 
algebraic notation for the parameter 7: 

(1) If c/ = 1 no construction is involved; we then write 7 = 0. 

(2) If c/ = 2 there are 3 distinct modes of contact of the two cubic cells, generated by face-wise, edge-wise, or 
vertex-wise joining of the cells; these modes specify the cluster geometry completely. We write 7 = </> (face-joining); 
or 7 = e (edge-joining); or 7 = w (vertex-joining) respectively. 

Since for instance face-joining <j) is not a geometrically unambiguously defined operation in the CA lattice space, we 
introduce more elementary construction symbols which specify unique and independent operations; the elementary 
(and in this specific case irreducible) independent operations are the joinings along the lattice axes: x axis, (f>x', y axis, 
(py; and z-axis, (p^- (We can join the second cell to the first cell in the positive or in the negative x-direction; however, 
the two resulting configurations are superposable under translations along the a;-axis; they count as 'geometrically 
equal'). 

These construction symbols can be combined by a first operation, of addition (+) 

(t>x + (l^y + (l^z = (t>; 

the sign (+) is read 'or': face-joining is face joining along the .x-axis, or along the y-axis, or along the z-axis, and this 
exhausts the possible alternatives. Denote the number of cell -faces in the x - direction by ^pxi etc; we then have ip^ (= 
<^j/ —'^z) — '/'/S (=2). The multiplicity of the geometric configuration generated by (j)x is 1 (=iy92:/2, the division by 2 
being due to the geometrical equivalence of the joining along the 'positive' or 'negative' face). The multiplicity of the 
geometric-cluster class G(2, 2; (f) is then the sum of the multiplicities of the component elementary geometric-cluster 
classes G{2, 2; (j)^), etc. This property of additivity of multiplicities is an instance of the following obvious property: 
If a general construction procedure 7 is the sum (+) of independent elementary construction procedures, 71, 72, 
. . ., then the multiplicity of the cluster-geometry is the sum of the multiplicities of the corresponding elementary 
cluster-geometries: 

7 = 71 + 72 + ••• + 7n : 9{a,c-n) = g{a,c;^i) + 3(0,0; 72) +... + 3(0,0; 7„) . (Al) 

A similar break-up of the operation of edge-joining e holds, e = Cx +ey +6^; where Cx is the more elementary (though 
not irreducible) operation of joining two cells along an edge parallel to the x-axis. We have, with obvious notations, 

(= Cj/ = £«) = e/3 (=4) different edges parallel to the x-axis, and, as in the case of face-joining, half of these edges 
(=6x72) produce geometrically distinct configurations (the joining along opposite edges with respect to the centre 
of the cube generates geometrically equal configurations; the joining along edges belonging to a same face produces 
geometrically distinct configurations). 

Similar considerations hold for vertex-joining. 

(3) If Cf = 3, first form a 2-cell cluster, joining two cells as under (2); the third cell is then attached to the 2-cell 
cluster. A sequence of two joining operations is indicated algebraically by a multiplication sign (.) between these 
operations. For instance, the combined operation ey.ex symbolizes the construction of 3-cell clusters, a pair of cells 
being joined along an x-edge; the attachment of the third cell is along a y-edge. Through addition and multiplication 
of elementary construction operations more complex construction schemes can be generated. 

For c/ = 3 we have 6 main combinations of joinings, {(j) . (p, e . e, v . v, (j) . e, e . v, v . (p); each main combination 
gives rise to directional variants which are described in terms of the elementary operations {4>x, etc). As an instance, 
take the face-wise joining of all 3 cells, the common faces being parallel; this operation is symbolised by 4>x ■ <f>x + 

(py . (py + . (Pz- 

In the construction of cluster geometries we have to take account very often of special clauses in the combination 
of the elementary operations. For instance, among the different edge-joinings of 3 cells, consider the joining along 

x-edges combined with the joining along y-edges, {cy . Cx), under the extra constraint that both edges have no 
vertex in common. Such extra clauses are symbolised by appropriate subscripts to the construction symbol. In our 
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illustration the absence of a common vertex is indicated by the subscript —v, so that the construction symbol becomes 
{^y^x)-v In a similar fashion the symbol {ey.ex)+v points out that the two edges are required to have a common 
vertex (subsript +v). We have 

(This notation is consistent with the symbols for the elementary operations; in cpx the subscript x indicates a constraint 
to the general face-joining operation (j)). 

The trivial cluster equivalence class C(l) contains a single cluster geometry, G(l, 1; 0) (single non-empty cell). The 
corresponding binding energy is zero, 

G(1,1;0) : E{1,1;0) = , g{l,l;(D) = 1. {A2) 

C. Cluster equivalence class C(2) 

The cluster equivalence class C(2) is (a) either a collection of c = 2 adjacent cells containing one particle each; (b) 
or it is a single cell, c =1, containing both particles. 

(a) Under the first alternative we have all three modes of joining the two cells, G{2, 2; j), j = (p, e, or v: 

(i) Common face: geometry G(2,2;0). 
Multiplicity : 

ff(2,2;<^) = 5(2, 2; + g{2,2;cl)y) + g{2,2; cj^y). 
Since g{2, 2; (p^) = ¥'x/2 (cf above), we have 

G(2,2;.^) : £(2,2;<^) = V^, , g{2,2-4>) = f = 3. (A3) 

(ii) Common edge: geometry G(2,2;e). 

Proceeding as under (i) (with substitution — > e, and — * e), we have g{2,2;ex) = ex/2, and hence 

G(2,2;e) : i?(2,2;e) = , g{2,2;e) = | = 6. {AA) 

(iii) Common vertex: geometry G(2,2;w). 

In (i) substitute (p ^ v, and (p ^ u, and consider vertices lying on the main diagonals of the cubic cell; take acount 
that the joining in opposite directions along a diagonal generates geometrically equal configurations. Therefore 

Gi2,2;v) : E{2,2;v) = V, , gi2,2;v) = ^ = 4. {A5) 

(b) The second alternative is trivial 

G(2,l;0) : E{2,1;<D) = V^, 5(2, 1;0) = 1- {A6) 

D. Cluster equivalence class C(3) 

This class C{3) gives rise to 3 broad categories of geometrical clusters, (a) G(3, 3; 7), (b) G(3, 2; 7), and (c) G(3, 1; 0). 

(a) Category G(3, 3; 7) contains a first variety of 3 geometries of type G(3, 3;j . j), obtained by attaching one cell 
to each of the 3 geometries G(2, 2; j), {j designates an operation (j), or e, or v, as under G(2); if the first joining is a 
face joining, then so is the second; etc); this variety of geometrical clusters is constrained to form three aligned cells. 

(i) Pair of common opposite faces: geometry G(3; 3; (0 ■(f))-e) (three cells aligned along a lattice axis); subscript — e 
indicates that the two faces of the cubic cell have no common edge; hence they are parallel. 

The construction procedure is explicited as follows (</> ■(p)-e = <Px ■ <Px + fpy ■ (t^y + ^t'z ■ 4'z', hence the multiplicity 
from Eq (Al): 



G(3,3;(.^. .^)_e) : E{3,3; {cl> . <f>)-e) = 2V^ , 5(3, 3; (,/. . 0)_e) = f = 3. 
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(ii) Pair of diametrically opposite eommon edges: geometry G(3, 3; (e . e)_0) (three cells aligned along diagonals of 
one face); subscript —0 indicates that the two edges of the cubic cell do not belong to a common face; this implies 
that they are symmetric with respect to the cell centre. 

G(3,3;(e . e)_^) : i;(3, 3; (e . e)_^) = 214, 5(3, 3; (e . e)_^) = | = 6. {A8) 

(iii) Pair of diametrically opposite common vertices: geometry G{3,3;{v . w)-^) (three cells aligned along cell 
diagonal); as under (ii) subscript —(p stresses that the two vertices do not belong to a common face; hence they are 
symmetric with respect to the cell centre. 

G(3,3;(« . «)_^) : E{3,3; [v . v)_^) = 2K , g{3,'S; {v . v).^) = ^ = 4. (A9) 

Category G(3, 3; 7) includes a second variety of 3 geometries of type G(3, 3; j . j), this time under the extra constraint 
of non-alignment of the three cells. (The two constraints, alignment, and non-alignment of the three cells, joined by 
a same operation j, manifestly exhaust all alternatives of type G(3, 3; j . j)). 

(i') Pair of non-parallel common faces: geometry G(3,3; (0 . (j>)+e) (L-shaped configuration); subscript +e points 
out that the two faces have a common edge. 

Start with a central cube, and add a pair of cubes having faces in common with the central cube; there are ip{ip—l)/2 
different pairs, among which (p/2 pairs belong to alternative (i) already listed. All of these configurations have same 
energy. 

2 

G(3,3;(<^. </))+e) : £(3, 3; (<^ . <^)+e) = 2V^ + V, , ^(3, 3; (<^ . <^)+e) = ^ - = 12. (^10) 

(ii') The case of non-diametrically opposite common edges gives rise to 3 geometry classes. 
(q) Parallel joining edges: The construction scheme is 

(e . e)+^ = (ca, . ex)+<p + (e^ . 6^)+^ + (e^ . 6^)+^, 

where the subscript +</> signifies that the two parallel edges belong to a same face of the cube. The multiplicity of 
{ex-^x)+(t> is (cf i') £3,^/2 — Cx, with = e/3. From Eq (Al) we obtain the total multiplicity. 

G(3,3;(e.e)+^) : £(3, 3; (e . e)+^) = 214 , ^(3, 3; (e . e)+^) = ^ - e = 12. (^11) 

(/?) The two joining edges are distinct, e, e' (mutually orthogonal), and have no common vertex (subscript —v); 
(alternatively we may say that the two edges do not belong to the same face). Hence the construction scheme is 

(C . 6 )— 1, (Ca; . Cy^ —y -\- (Cy . € z) —y "t" (Cjg • ^x)—V 

The multiplicity of {cx ■ ey)-v is (cf i') €xey/2. Hence 

G(3,3;(e.e')-„) : £(3, 3; (e . e')-.) = 2V, , ^(3, 3; (e . e')-.) = ^ = 24. (^12) 

(7) The two joining edges are distinct (mutually orthogonal), and do have a common vertex (subscript +v) (or they 
do belong to a same face). Under this alternative each cell shares an edge with each of its two neighbour cells. The 
construction scheme is 

(e . e )-|-u = ' ^y)+u ~t~ (^y • ^z^+v ~i~ (^z * ^x)+u* 

The multiplicity of the geometry generated by (e^ . 6^)+^ is as under (/?), namely ex€y/2. 

G(3,3;(e . e')+.) : £;(3, 3; (e . e')+.) = 2^ + ^(3, 3; (e . e')+.) = - = 24. {A13) 

(iii') Pair of non-diametrically opposite common vertices. These configurations fall into 2 geometry classes. 

(a) Two vertices lying on same edge: geometry G(3, 3; {v . v)+e) (subscript +e indicating that the vertices are 
joined by a common edge). We decompose this construction scheme as follows 

{V . v)+e = {v . z;)+e, + {v ■ v)+e^ + {v . v)+e,, 
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with the notation {v . v)+ex indicating that the two vertices are constrained to lie on an edge, etc. The multiplicity 

of {v . w)+e» is ex = 3. 

G{3,3;{v .v)+e) : E{3,3; {v . v)+e) = 2V, , g{3,3; {v . v)+e) = e = 12. (AU) 

Two vertices belonging to different edges of same face (vertices belonging to different faces are diametrically 
opposite and have been dealt with under (iii) A9). It is convenient to describe this alternatively by saying that these 
vertices belong to the same diagonal J of a face (constraint indicated by subscript +6): G{3, 3; {v . v)+s)- If we denote 
by (5j. the collection of the diagonals of the faces normal to the x-lattice direction, etc, the relevant construction 
scheme becomes 

{v . v)+s = {v . v)+s^ + {v . v)+s^ + {v . v)+s^. 
Since the multiplicity of {v . v)+Sx is readily seen to be 2(fix = 4, we have 

G{3,3;{v .v)+s) : E{3,3; {v . v)+s) = 2K , g{3,3; {v . v)+s) = 2<f = 12. (A15) 

Finally, category G(3, 3;7) generates a third variety of 3 geometries of type (7(3, 3; j . j'), where j and j' are of 

different nature (if j stands for face-joining, j' must be either edge-joining or vertex-joining, etc). 

(iv) Face and edge-joining: geometries 0(3, 3; e . (j)). These configurations fall into 2 geometry classes. 

(a) Edge joining on 'small side' of the face-joined box: G{3, 3; e\\^ . (p) (where ey^ indicates that the edge is parallel 
to the face <j)). Construction scheme: 

e||^ . (j) = By . +ez . (j)x +ez . (j)y + . (j)y + ■ (l)z +ey . (j)^ ■ 

The multiplicity of the geometries resulting from Cy . (j)x is ipx^x/^ (=4), hence 

G(3,3;e||^ . </>) : iJ(3, 3; e,,^ . </>) = V^, + V, , <?(3, 3; ey^ . </>) = ^ = 24. (A16) 

(/3) Edge joining on 'long side' of the face-joined box: G(3, 3; e_i_^ . <p) (where e±^ indicates that the edge is normal 
to the face (p). Construction scheme: 

e_L0 . (j) = . +ey . (j)y + . . 
Since the number of edges of the 'long side' is 2ex, the multiplicity of ^x^x is Vx^x (=4). 

G{3,3;e±^ . (j)) : Ei3,3;e±^ . (j)) = + Ve + Vy , ^(3, 3; . = ^ = 24. (A17) 

(v) Face and vertex-joining (third cell attached to an outer free vertex of box): geometries G(3, 3;v . (j)). Construc- 
tion scheme 

V . (j) = V . (j)x + V . 4>y + V . (j)z- 

The multiplicity oi v . (px is !^<Px/2 

G{3,3;v . (j)) : ^(3,3;^;..^) = V^, + Vy , g{3,3;v . cfy = ^ = 24. (A18) 

(vi) Edge and vertex-joining (two cells joined vertex-wise; third cell attached to an edge not incident to the joining 
vertex): geometries G(3, 3;'y . e). Construction scheme: Attach second cell to a vertex of the first cell {y/2 distinct 
alternatives, according to G(2, 2; w)); the third cell is to be added such as to have a common vertex (but not a common 
face) with the two cells; we have 3/4e (=9) alternatives on each of the two cells. Hence the multiplicity 

G(3,3;w.e) : E{3,3;v . e) = Ve + Vy , g{3,3;v . e) = ^= 72. (A19) 

(b) Category G(3, 2; 7) generates essentially the variety of geometries G(2, 2; 7), with one difference. In the G(2, 2; 7) 

case the two cells arc equivalent and indistinguishable; in the G(3,2;7) case, one cell contains 2 particles, while the 
other contains only one particle, so that the translational symmetries do not hold. The multiplicities are therefore 
twice the multiplicities encountered under G(2,2;7). 
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(i) Common face: geometry G{3,2;<p). 

G{3,2;<P) : E{3,2;<j,) = 2V^ + g{3,2;<j,) = ^ = 6. {A20) 

(ii) Common edge: geometry G(3, 2;e). 

G(3,2;e) : E{3,2;e) = 2Ve + Vo , g{3,2;e) = e = 12. {A21) 

(iii) Common vertex: geometry G{3,2;v). 

In (i) substitute (j) v, and consider vertices lying on the main diagonals; take acount that the joining in opposite 
directions along a diagonal produces an equivalent configuration under combinations of translations. Therefore 

G{3,2;v) : E{3,2;v) = 2Vy + ^(3,2;^;) = v = d>. {A22) 

(c) The final alternative of all 3 particles in the same cell generates the trivial geometry 

G(3,l;0) : ^(3,1;0) = 3^ + 3Ay, g{3,l;$) = 1. {A23) 

This completes our list of distinct geometrical configurations associated with the cluster classes C(l), C(2) and 
C(3). 

The systematic construction method of the configurations extends naturally to clusters of arbitrary size a. However, 
with increasing a the number of alternatives increases exponentially, so that a detailed enumeration becomes rapidly 
prohibitive. This suggests that an asymptotic approach, with a taken as the 'large' parameter, should be substituted 
to the detailed enumeration procedure of this Appendix. Common experience with asymptotic expansions suggests 
that results precise enough for the statistical purposes of Section 5 should become available even for fairly low a 
values. 
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FIG. 1: Initial spatial configuration (t=0) for head-on collision of two identical nuclei, T, I, of mass 150. Direction of 
propagation: a;-axis. T and I separated by empty layer of cells x = Q. 
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FIG. 2: Configuration of firagments at timesteps t — 20 (panel a), 40 (b), 50 (c), and 60 (d). Reconstructed symmetrized 
distribution of nucleons for timestep 40 (b'). 
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[(A(t)-A)/(A(0)-A)]'' 
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ln(A (t)-A^rn(A (0)-A ) 
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FIG. 3: Time-dependence of size of central cluster: [{Ac{t) - Ao)/{Ac{0) - Ao)]'' , k = 0.5, 0.4, 0.3, 0.2, 0.1 (thin curves), 
and {hi{Ac{t) — Ao)/{\n.{Ac{0) — Ag) against time (heavy curve). Impact energy: 3.75 MeV per nucleon). 
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FIG. 4: Cluster distribution log A'^ - log a derived from simulated 4n detectors, (a) Eimp/A = 0.307 MeV per nucleon: r = 
2.56; (b) Ei,np/A = 3.973 MeV per nucleon: r = 2.64. The numbers of clusters are normalized to the numbers of runs. 
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ENERGY/NIUCLEON (MeV) 

FIG. 5: Formal CA caloric curves, T-Eimp/A, for different definitions of temperature: (a) T'„(c)) nucleon-gas in compound 
nucleus; (a') Tn((7)i nucleon-gas in CA universe; (b) Tc(c), cluster-gas in compound nucleus; (&') T£.((7) ;cluster-gas in CA 
universe; (c) ri3/22, reaction equilibrium C(l) + C(3) = 2C(2); (c') T3/12, reaction equilibrium C(3) = C(l) + C(2); (c") 
^2/11, reaction equilibrium C(2) = 2C(1). Individual points: experimental estimates from Au-Au fragments (d) NuPECC 
interpretation of data (full squares); Trautmann interpretation (e, full dots), (f, crossed squares), (g, double crosses). 



